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Abstract. The two- layer disk models of Chiang & Goldreich (1997, henceforth CG) and its derivatives are popular 
among astronomers because of their simplicity and the clear predictions they make for the SEDs of T Tauri 
stars and Herbig Ae/Be stars. Moreover, they can be computed quickly, which is a great advantage when fitting 
observations using automated procedures. In this paper we wish to assess the accuracy and reliability of 2-layer 
models, by comparing them to detailed vertical structure models with accurate I+ID radiative transfer. We focus 
on the shape of the SED, and the predicted height and "flaring index" of the disk. We flrst consider models 
where scattering is set to zero. We find that 2-layer models overestimate significantly the near-infrared flux, and 
we suggest a simple way of correcting this effect, at least in part. At longer wavelengths, the SED of two- layer 
models often show a two-bump structure, which is absent in 1-l-lD models. Nevertheless, overall agreement is 
reasonably good, and the differences are in most cases within 30%. At (sub)-mm wavelengths the differences may 
even be less. The shape of the disk, as measured by its pressure and surface scale height and by the fiaring angle 
are also well reproduced by two-layer models. When scattering is included in the 1-l-lD models, the differences 
become larger, especially in the near-infrared. We suggest simple ways to include scattering in two-layer models 
and discuss their reliability. We do not compare the two-layer models to full 2-D/3-D models, so the conclusions 
remain valid only within the annulus-by-annulus approximation. 
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1. Introduction 

Models of disks around pre-main-sequence stars are cur- 
rently used to make predictions of a number of observ- 
able quantities and to compare them to observations. The 
complexity of such models varies from so-called power- 
law disks, where all the relevant disk properties are de- 
scribed by power-law of radius and changed independently 
(see, for example, Beckwith et al.^990j), to models where 
the thermal and geometrical properties are computed self- 
consistently, under th e assu mption of hyd rostatic eq uilib- 
rium (D'Alessio et al. Il998t Malbe t et a l. Il99ll200li Bell 
et al. ll997lll999t Dullemond et al. l2002|) . A large effort is 
currently being made in improving the treatment of radia- 
tion transfer, by developing 2D codes that deal efficiently 
with the very large optical depths that characterize such 
disks (e.g. Bjorkman & Wood JOQl, Dullemond .20021 . At 
the same time, it is clear that the need of simple mod- 
els, that can be easily used in interpreting the observa- 
tions, remains. This explain the success of the two-layers 
disk schemat ization proposed some years ago by Chiang 
& Goldreich lll997t henceforth CG97). 



Because the CG97 model and its derivatives 
(e.g. Chiang et al. l200l[^ Dullemond, Dominik & Natta 
.2001; Lachaume et al. l2003|) are often used in interpret- 
ing disk observations, it is important to estimate their 
reliability and accuracy by comparing them to more de- 
tailed and self-consistent model calculations. At present, 
self-consistent 2-D and 3-D models are still too cumber- 
some for such a comparison. However, there are a number 
of ID vertical structure models based on detailed radia- 
tive transfer, which can be used as templates to assess the 
reliability of 2-layer models. 

In this pape r, we will make use of the model of 
Dullemond et al. l2002l Chenceforth DZN02). In this model 
the vertical temperature profile was solved using a de- 
tailed plane-parallel 1-D vertical radiative transfer proce- 
dure, which treats properly the slanted penetration of the 
stellar radiation (often referred to as 1-1- ID approxima- 
tion). Using the resulting temperature profile, the density 
structure is re-computed after each iteration step by inte- 
gration of the equation of vertical hydrostatic equilibrium. 
In a follow-up paper (Dullemond & Natta 200^ scattering 
of the stellar light by the dust particles was included in 
the model. For the case of zero albedo this model is in fact 
identical to the DZN02 model, and we will use DN03 as 
our "standard 1-1- ID disk structure model" against which 
we will compare the two-layer models. 
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The goal of this paper is three-fold. First, and most 
importantly, we compare an improved version of the CG97 
model to the the 1+lD vertical structure model of DN03. 
This initial comparison is done for dust opacities in which 
scattering is switched off. Then we investigate how big 
the effects of scattering are in comparison to the typical 
errors of the improved CG model. Finally we discuss the 
possibility of using a simple recipe to include scattering 
into the CG model. 



2. Description of the models 

In this section we give a short description of the improved 
CG97 model we adopt, and of the 1-1- ID vertical struc- 
ture model. Both families of models neglect viscosity, and 
treat the radial dependence of the surface density as a free 
parameter. Consequently, they also do not include viscous 
heating which, depending on the value of the viscosity pa- 
rameter Of, could be of importance in the inner parts of 
the disk. 

2.1. Improved CG97 model ("CG+") 

The central element of the two-layer Chiang and Goldreich 
model is that, due to the very shallow incident angle 
a ^ 1, the impinging stellar radiation is absorbed en- 
tirely in the very tenuous upper layers of the disk, high 
above the disk's photosphere. While this layer has an op- 
tical depth of unity for stellar radiation at that grazing 
angle (typically around a = 0.05), it is very optically thin 
in the vertical direction (rverticai — ct). Energy conserva- 
tion requires that this absorbed radiation is re-emitted as 
infrared radiation, half of which is emitted away from the 
disk, while the other half is emitted towards the disk mid- 
plane. This downwards emitted flux is absorbed by the 
disk interior and re-emitted once more in the infrared, 
though this time at longer wavelengths. This typically 
gives two components in the SED: a warm optically thin 
component with dust emission features in the near- to mid 
infrared and a cool thermal blackbody component from 
the disk interior at longer wavelengths. If the disk is opti- 
cally thick to the radiation from the midplane, then both 
components are equally strong, as the cool component is 
a reprocessed version of the downward radiated surface 
emission. The temperature of the surface layer is fixed: 
it is the optically thin dust temperature. The emissivity 
of the surface layer is determined by the surface density of 
the surface layer, which is calculated from energy conser- 
vation: the total emission should equal the total absorbed 
flux. The emission from the interior is determined by the 
temperature of the disk interior and the disk's surface den- 
sity. In the simple case of a grey opacity and high optical 
depth this so-called midplane temperature is determined 
by equating the thermally emitted flux F = aT^ to the 
downwards directed surface flux. In general, Tj is com- 
puted using mean opacities both for the heating and for 
the cooling radiation. 



The original paper of Chiang & Goldreich (Il997l) de- 
scribes the two-layer flaring disk model in stages. At first, 
a grey opacity constant surface density model is presented 
for which the flaring index (ilog(i/s)/dlog(i?) = 9/7, 
where Hg is the height of the surface. Then a set of im- 
provements are listed in order to account for effects of 
low optical depth at vari ous w avelength regimes. In a 
later paper (Chiang et al. l200lj) the model was refined 
to include more realistic opacities and, more importantly, 
the self-consistent (and numerically stable) determination 
of the fiaring angle, which is crucial f or ene rgy conser- 
vation. Dullemond, Dominik & Natta l)200ll henceforth 
DDNOl) subsequently added the emission and structure of 
the disk's inner rim, and included effects of self-irradiation 
and self-shadowing. The complete set of model equations 
is listed in DDNOl. In the rest of this paper this improved 
Chiang & Goldreich model will be referred to as the CG"*" 
model. We will not address here issues related to the inner 
rim. 

2.2. Vertical structure models 

The reference model against which we compare the CG"*" 
model is the 1-l-lD vertical structure model of DZN02 and 
its improved version DN03. The radiative transfer is done 
in a two-stage procedure: first the primary stellar photons 
are followed as they penetrate the surface layers of the 
disk and get absorbed. Then the re-emitted infrared radi- 
ation is solved using a 1-D vertical frequency and angle- 
dependent radiative transfer code applied at each radial 
annulus. In a third stage the vertical pressure balance is 
solved, under the assumption that the gas temperature 
equals the dust temperature and that the gas-to-dust ra- 
tio is constant. This three-stage procedure is then iterated 
until a converged solution is reached for the complete tem- 
perature and density structure of the disk. During each it- 
eration, the flaring index has to be recomputed in order to 
compute the fiaring angle ( "grazing angle" in the language 
of Chiang & Goldreich) self-consistently. This is done by 
estimating the surface height of each annulus by finding 
the height Hg above the midplane where the direct stellar 
radiation is extincted by exp(— 1). The double-logarithmic 
derivative of this with cylindric radius R is then the 
fiaring index from which the flaring angle can be computed 
(see DZN02). 

Once the iteration procedure has converged, we have 
obtained a solution for the dust temperature Tdust (R, Z) 
and the disk's gas-|-dust density p{R, Z). 

DN03 included scattering in the following way. Stage 
1 of the DZN02 model (the irradiation of the disk by pri- 
mary stellar photons) was replaced with a 1-D Monte- 
Carlo code. Primary stellar photons are followed as they 
penetrate into the surface of the disk and scatter around 
until they either get absorbed, or escape to infinity. If a 
photon is absorbed, it leaves behind its energy, which is 
then used as a source term for the radiative transfer of 
the second stage. In the second stage, which takes care of 
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the re-emitted infrared radiation from the disk, scattering 
is not included. This is a valid approximation if the dust 
grains are small enough (typically not larger than 0.5 /im). 

3. Comparing models without scattering 

As a first step we compare the CG^ model and the l-flD 
vertical structure for a test case without scattering. In this 
way the physical assumptions of both models are the same, 
so that in principle we expect both models to produce very 
similar SEDs. The comparison is first done for a single 
annulus at various radii and optical depths, irradiated by a 
T Tauri star or by a Herbig Ae/Be star. We then compare 
the SEDs of models of entire disks, as well as intensity 
profiles at some selected wavelengths. 

3.1. Single annulus 

As a first example we consider an annulus of a disk at 
1 AU from the central star, with width 0.01 AU (i.e. an 
annulus between 1.00 and 1.01 AU). We take the central 
star to be a Herbig Ae/Be star (Af, = 2Mq, = 2Rq, 
Tcs ~ lOOOOK). To simplify matters, we assume the spec- 
trum of the star to be a blackbody spectrum with T = Tcs. 
The grazing angle (i.e. the angle at which the stellar ra- 
diation enters the disk's atmosphere) is computed self- 
consistently in both CG+ and 1+lD models, but for the 
purpose of the comparison described in this section we as- 
sume it to be fixed to a = 0.05. We also assume, for these 
annulus tests, that half of the surface of the star as seen 
from the top of the annulus is covered by the very inner 
parts of the disk (see CG97). This effectively reduces the 
stellar flux by a factor 0.5. The absorption cross sectio n 
is that of a silicate grain of 0.1 fim (Draine & LeeTosi), 
with zero albedo. Different choices of the opacity will not 
change the results of this paper in any significant way, 
unless explicitely noted. The surface density of the disk is 
such that the total vertical absorption optical depth (from 
z = —oo to z — oo) at 550 nm is T^QOnm ~ 10. The mean 
molecular weight is assumed to be /i = 2.3. 

In Fig.^the SED of this annulus is shown for both the 
vertical structure model and the CG"'' model. The differ- 
ence between the models is less than 30% at wavelengths 
A > 8/im. At shorter wavelengths the GG+ model has 
a bump that is not present in the SED of the vertical 
structure model. The reason can be traced back to the 
surface temperature assumed in the CG^ model, which 
is that of grains in an optically thin medium. In reality, 
however, only a fraction of the grains in the surface layer 
see the unextincted stellar light. More than 36 % of the 
material in the surface layer lies at optical depths larger 
than unity (along the grazing incident path of the stel- 
lar radiation) . It is hard to assign a single temperature to 
the material in the surface layer. But a reasonable aver- 
age temperature could be the temperature at an optical 
depth of T = 0.7, which is the location where exactly half 
of the stellar radiation has been absorbed. The dashed 
line in Fig. ^ shows the SED as computed by the CG+ 
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Fig. 1. The SED for the single-annulus test case nr. Al 
(Herbig Ae star) at 1 AU, compared to the GG^ model in 
its two proposed variants. 
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Fig. 2. The temperature structure for the single-annulus 
test case Al (Herbig Ae star), compared to the GG"*" model 
in its two proposed variants. The transition from the sur- 
face to the midplane temperature in the GG"*" models oc- 
curs at Hg, i.e., where the optical depth along a ray at a 
grazing angle a = 0.05 at stellar wavelengths is r, ~ 1. 
The tick mark on the solid curve shows the location of Hg 
in the 1-l-lD case. 



model in which the surface temperature is evaluated at 
an extinction of exp(— 0.7) — 0.5. 

In Fig.[21the temperature structure of the same annu- 
lus is shown as a function of the vertical coordinate Z/R, 
where Z is the height above the midplane, and R the ra- 
dius from the star (in this case R = lAU). High above the 
surface of the disk the temperature is constant, and vir- 
tually equal to the optically thin dust grain temperature 
at that distance from the star. The slight difference is due 
to the fact that in the 1-D model the dust high above the 
disk is heated not only by the direct stellar radiation but 
also by the thermal emission of the disk below it. 
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Ty 
10 
10 
300 



Al 2 2 10000 1. 
A2 0.5 2 4000 0.1 
A3 2 2 10000 50. 

Table 1. The parameters of the single- annulus models 
used in the comparison. The first column is the identi- 
fication used in this paper. Columns 2-4 are the stellar 
parameters (in units of Mq, Rq and K) and columns 5- 
6 are the annulus parameters. The radius is in AU. All 
model annuli have a width AR/R — 0.01 and a grazing 
incident angle of stellar radiation of a = 0.05. 



As one goes towards smaller Z the temperature starts 
to decrease due to extinction of the direct stellar flux 
by the material of the disk. This is the location of the 
disk's "hot surface layer" which produces the dust emis- 
sion features. The small tick mark on the temperature 
curve denotes the surface height Hs where the grazing 
optical depth is unity, i.e. where most of the direct stel- 
lar radiation is absorbed. Continuing downwards from the 
surface layers, the temperature slope levels off again be- 
cause the disk's own infrared radiation becomes competi- 
tive with the direct stellar radiation in the thermal balance 
equation of the dust grains. This is where the disk's pho- 
tosphere is located and where most of the thermal contin- 
uum emission is produced (the midplane emission, in the 
terminology of CG models). 

Over-plotted over the smooth temperature curve is 
the two-layer CG"*" model, in its two variants mentioned 
above. The sharp step in the temperature profile is where 
the surface layer of the CG+ model starts: z = Hs- The 
main difference in the two variants of the model is the tem- 
perature of this surface layer. For one variant, the surface 
layer temperature equals that of optically thin dust grains. 
For the corrected variant this temperature is lower, and is 
closer to the temperature of the vertical structure model 
at the location of the tick mark. 

Note that Hg is different in CG+ and 1+lD models. 
Most of difference is due to the fact that in the 1-1- ID 
models the density in the disk surface deviates from the 
Gaussian shape predicted by isothermal models (see, for 
example, Fig. 3 of DN03). An additional, smaller differ- 
ence is due to the fact that CG"*" models use the mean 
opacity averaged over the stellar flux to define Hg, while 
the 1-l-lD models use the wavelength-dependent opacity. 

Tabled lists the parameters of two other annuli. Test 
A2 is the T Tauri analog of Al, with lower stellar mass and 
temperature, and the annulus taken at a smaller radius (to 
obtain similar temperatures in the disk). In Fig.|3]the SED 
of this test case is shown. The results are qualitatively 
similar to the case of the Herbig Ae star (test Al). In 
this particular example, however, the exp(— 0.7) correction 
of the surface temperature calculation in the CG"*" model 
does not make a major improvement, although also in this 
case the overall SED of the 1+lD model is slightly better 
reproduced when this factor is taken into account. 
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Fig. 3. The SED for the single-annulus test case nr. Al 
(Herbig Ae star) at 1 AU, compared to the CG^ model in 
its two proposed variants. 

As shown above, in many cases the CG'^ model is rea- 
sonably well in agreement with the l-flD vertical struc- 
ture model. Yet, there are cases in which the differences 
are appreciable. We show an example in Fig. 0] which plots 
the SED of the A3 annulus. This annulus has Ty — 300, 
and is at 50 AU from our Herbig Ae star. It is clearly seen 
that in this case the difference between CG+ and 1+lD 
predictions is rather large, though still within a factor of 2. 
The explanation for the breakdown of the CG+ model in 
this case is likely the use of mean opacities to compute the 
midplane temperature. As noted by DZN02, the midplane 
of the disk can cool down by leaking radiation at long 
wavelengths (where the optical depth of the disk becomes 
smaller than unity) , more than expected on the basis of a 
mean-opacity analysis. This is clearly seen in Fig.[Sl which 
shows the temperature structure of the A3 annulus. The 
1-l-lD temperature profile, which goes from values lower 
than predicted by CG^ in the midplane to larger ones at 
intermediate z results in lower fluxes around 100 /im and 
higher fluxes in the mid-infrared, so that the two-bump 
shape of CG+ practically disappears in the 1-l-lD model. 

3.2. Full disk models 

In this section we compare the SED and the physical 
structure of entire disks computed with the I+ID vertical 
structure model and CG+ models. We present three illus- 
trative cases: a disk around a Herbig Ae star, a T Tauri 
star and a Brown Dwarf. Table |2] lists the parameters of 
the three disk models. 

We assume that the disk is in hydrostatic equilibrium 
and compute the flaring angle is self-consistently. This is 
done by evaluating the flaring index 

d\g{Hg/R) 



(1) 



d\gR ' 

where Hg is the height of the Tgrazing = 1 surface of the 
disk. From the flaring index, the flaring angle can be com- 



DuUemond & Natta: An analysis of two-layer models for circumstellar disks 



5 



10 



10 



10 



10 



10" 



10 



6« 50 ■ 

i ■ 



<D -50 ■ 



R = 50 AU 
AR = 0.5 AU 
a = 0.05 
i = 



10 



Herbig Ae (A3) 




1 + 1D 
CG + 

CG + , exp(-0.7) 



100 



100.000 
10.000 
1.000 
0.100 
0.010 
0.001 



40 
20 


-20 



1000 



Full disk 


Herbig Ae (F1) 


(no disk rim) 




j 






\ 


' r 


1+1D \ 




CG+ \ - 


■ r 

- ■ A 

-1 — 1 1 1 / 1 1 1 — J — 1— ^ 


CG + , exp(-0.7) \ 



















10 



100 



1000 



A [/in 



Fig. 4. The SED for the single-annulus test case nr. A3 Fig. 6. The SED for the fuU disk model Fl (Herbig Ae 
(Herbig Ae star, high optical depth) at 50 AU, compared star), compared to the CG^ model in its two proposed 
to the CG"*" model in its two proposed variants. variants. 
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The temperature structure form model A3 

i?in i?out E(1A(7) " 



M. 



Mdis 



Fl 2 2 10000 1 300 1000 0.2 
F2 0.5 2 4000 0.1 300 100 0.02 
F3 0.1 1.3 2600 0.033 30 100 0.002 

Table 2. The parameters of the full disk models. The first 
column is the identification used in this paper. Columns 
2-4 are the stellar parameters (in units of Mq, Rq and K) 
and columns 5-6 are the annulus parameters: outer radius 
in AU, and the surface density E at 1 AU. All models have 
a power law index for the surface density of -1 (S cx 1/i?)- 
Column 8 is the mass of the disk in units of Mq as derived 
from the other parameters. 



puted directly: 



(2) 



For numerical stability reasons this fiaring index is always 
evaluated 2 grid points inwards of the point where it is 
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Fig. 7. Same as Fig. but now for model F2 (a T Tauri 
star). 



used (see Chiang et al. l200ll) . In these models we assume 
that the stellar fiux is not reduced due to occultation of 
the central star by the inner disk, as was assumed in the 
single-annulus models. 

Fig. shows the SED of model Fl. Over-plotted are 
the two proposed CG"*" models (the one with optically thin 
surface temperature and the one with the correction factor 
exp(— 0.7)). One can see that overall the models agree, 
but that the CG+ model clearly shows a double-bump 
structure between 20 /im and 100 fim which is not seen 
in the full vertical structure model. In Fig. |7| the same 
is plotted, but now for a T Tauri star (model F2). The 
disk is taken to be less massive than the Herbig star disk. 
The CG"*" model works really well for this case, with errors 
mostly within 10%. Finally, in Fig. |Slthe SED of a Brown 
Dwarf disk model is show (model F3). This disk is taken 
to have a much smaller outer radius and mass than the T 
Tauri case. Nevertheless, the disk remains very optically 



6 



DuUemond & Natta: An analysis of two-layer models for circumstellar disks 



10 



10 



-2 



Brown Dworf (F5) 



"s 10" 



15 10 



10 



10" 





1+1D 
CG + 

CG + , exp(-0.7) 



40 ■ 



20 - 
'-' 
20 - 

-40 ^ 

1 



0.50 



10 



100 



1000 



Fig. 8. Same as Fig. El but now for model F3 (a Brown 
Dwarf). 



thick, as in the other models. Again, errors are reasonably 
small, though slightly larger than for model F2. 

In Fig. El the height of the disk of model Fl (Herbig 
Ae star) and F2 (the T Tauri case) is shown as a function 
of radius, for both the I+ID model and the CG+ model. 
The pressure scale height is derived from the temperature 
of the disk at the equatorial plane. At very small radii Hp 
is virtually the same in CG"'" and l+lD models. At larger 
radii, Hp in the l+lD model drops below that of the CG"*" 
model, because the equatorial plane temperature of the 
l+lD model tends to be lower, as discussed in Section 
18.11 The difference, however, is always very small. 

The surface scale height in Hg tends to be larger 
in l+lD models, reflecting in part the deviation of the 
very upper layers from a Gaussian profile, as discussed of 
Section rm The degree by which CG+ models can repro- 
duce the results of 1-1-1 D ones depends on the dust opacity 
in combination with other model parameters, such as the 
disk mass. Here we just note that the difference is larger 
for the HAe star, and for larger radii. However, these dif- 
ferences remain small and we can conclude that, in all 
cases we have considered, the disk shape can be derived 
with good accuracy using the simple approximations of 
the CG"^ models. 

Fig. 1 101 shows the flaring index as function of radius for 
model F2. This quantity gives information on the amount 
of radiation emitted by the disk at any given radius, be- 
cause for R ^ the emission is proportional to ^. 
Although the run of ^ with R is similar in the two models, 
we note that ^ is a smooth function of radius in the CG"*" 
models, but not in the 1+lD one. In particular, the kink 
in the flaring index 0.3 AU is not due to numerical inaccu- 
racies but rather to the behavior of the midplane optical 
depth , which in this disk model drops below unity at long 
wavelengths at about that radius. The consequent cooling 
of the midplane below optically thick temperatures, dis- 
cussed in §3.1, is seen as a slight kink in Hg in Fig. El This 
effect, which depends on the disk and stellar parameters 




R [AU] 



0.50 



0.40 



0.30 



0.20 



0.10 



0.00 




100.0 



Fig. 9. Top panel: Pressure scale height Hp and surface 
height Hs as function of R for model Fl (Herbig Ae star). 
The solid line shows the results of the vertical structure 
model, the dashed line those of CG"*" ones. Bottom Panel: 
same for model F2 (T Tauri star). 
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Fig. 10. The flaring index of the disk ^ is plotted as a 
function of radius for model F2 (T Tauri star). The ver- 
tical structure model is shown as a solid line, while the 
CG+ is the dashed line. 
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Fig. 11. The brightness profile of model F2 (T Tauri star) 
as a function of radius at two different wavelengths. The 
I+ID model results are shown as a solid line, the CG+ 
model results as a dashed line. 
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Fig. 12. Example of the effect of scattering on the SED 
of a T Tauri disk (model F3). Only the thermal emission 
of the disk is shown in this figure. Scattered light would 
appear as a strong rise towards short wavelengths. 



and on the exact wavelength dependence of the dust opac- 
ity, can have important consequences on the stability of 
1-f-lD models. 

Fig. 111! shows the brightness profiles of the two model 
approaches at two wavelengths, 10 /zm and 1.3mm, where 
current instrumentation has the spatial resolution re- 
quired to resolve the disks. The differences between the 
models are very small, probably within the observational 
accuracy of most experiments. The 10 /im profile has a 
clear exponential cut-off at large radii due to the fact that 
the dust becomes too cool to emit strongly at mid-infrared 
wavelengths. The CG+ and 1-l-lD models produce virtu- 
ally identical emission at this wavelength, as seen also in 
the SED (Fig. EJ. At 1.3 mm the differences are slightly 
larger (up to 30%) , which reflects the fact that sometimes 
the 1-l-lD models can be slightly cooler in the disk's mid- 
plane than the CG"*" models, as discussed above. 

4. Models with scattering 

4.1. Effect of scattering on the disk 

The effects of scattering of the stellar radiation on the 
structure and SED of circumstellar disks was discussed in 
detail in DN03. The main effect is that the thermal in- 
frared flux from the disk is reduced at all wavelengths, 
with a larger effect at shorter than at longer wavelengths. 
This reduction takes place because the scattering reflects 
away from the disk part of the stellar radiation, which 
therefore cannot be reprocessed into infrared radiation. 
The reflected fraction of the incident flux can be evalu- 
ated analytically using the H-functions of Chandrasekhar, 
which were discussed in this context in DN03. For a mod- 
erate value of the albedo {oj ~ 0.5) and isotropic scat- 
tering the reflected fraction is about 25% (of order lj/2), 
and the reduction factor of the disk's thermal emission is 
about 0.75, when averaged over the infrared spectrum. For 



higher albedo the reduction is larger, while it decreases if 
the scattering is more forward-peaked. 

The reason why the reduction is stronger at short 
wavelength than at long wavelengths can be traced back 
to the fact the surface layers of the disk are affected more 
strongly by scattering than the interior. A photon that 
gets lost in the disk's interior may scatter many times and 
still get absorbed rather than escape. Therefore the deep 
interior layers of the disk will have a temperature only 
slightly below the temperature of a disk without scatter- 
ing. 

The SEDs of a T Tauri disk (model F2) with increas- 
ing albedo, computed as in DN03, are shown in Fig. ^1 
The effect of reddening of the SED as oj increases is 
clearly seen. This suggests that CG+ model-predicted 
SEDs, which tend in general to overestimate the amount of 
flux at shorter wavelengths with respect to 1-l-lD models, 
may be significantly in error when realistic dust models 
are considered. 

4.2. A simple recipe for scattering in CG~^ models: is it 
possible? 

Although the qualitative effect of scattering can be eas- 
ily understood "a priori" , it is more difficult to estimate 
quantitatively how the infrared spectrum of the disk is 
affected at each individual wavelength. Still, it is impor- 
tant to explore the possibility of including scattering in 
the CG"*" models in some simple way, and to assess the 
reliability of any such model by comparing its predictions 
to those of the more accurate 1-l-lD models. 

A first recipe is based on the following considerations. 
In CG+ models, the stellar radiation is intercepted by the 
disk surface, which will re-emit (part as scattered light 
and part as reprocessed radiation) 1/2 of it toward the 
observer and 1/2 toward the midplane. This entire second 
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1/2 will be absorbed by the midplane, whose temperature, 
therefore, will not depend on u. Of the 1/2 emitted toward 
the observer, a fraction lj is in the form of scattered light, 
and consequently only a fraction (l-oj) is in the form of 
thermal reprocessed flux at infrared wavelength. Since the 
temperature of this surface layer dust depends only on its 
absorption cross section, one can compute the SED by 
combining the emission of the midplane of a CG^ model 
with Lo — Q with the surface emission reduced by a factor 
(l-w/2). This recipe neglects multiple scattering and any 
wavelength-dependent effect. However, for low albedo it 
reproduces well the overall reduction of the infrared ex- 
cess, and it is based on very simple arguments. 

A second possibility, suggested also by CG97, is to ap- 
ply the classical result that the fraction of incident photons 
absorbed by a plane-parallel slab is reduced by a factor 
rj = ^1 — Lo. The midplane will then have a temperature 
(in the optically thick case) lower than T{u) ~ 0) by a 
factor ?7^/'* = (1 — ui)^^^ . The surface grain temperature is 
again unchanged, but the infrared emission of the surface 
needs to be reduced by a factor rj. For not too high albedo, 
also this recipe will reproduce the overall reduction of the 
infrared excess found by DN03. 

We show a comparison of these two prescriptions with 
the correct results of DN03 in Fig.^1 The top panel shows 
the comparison for the Herbig Ae star disk (Fl), the bot- 
tom panel for the T Tauri star disk (F2). Each curve shows 
the fractional difference of the SED (reprocessed radiation 
only) of CG+ models, as labeled, to that of the comparison 
model, which is the l-flD one with oj — 0.5 (isotropic scat- 
tering). The solid line shows the results for CG"*" models 
where scattering is ignored. One can see that these mod- 
els overestimate the flux at all wavelengths, and that the 
discrepancy can be larger than a factor of two (for the 
Herbig Ae disk) between about 2 and 5 /im. At longer 
wavelengths the agreement improves, and becomes quite 
good below about 100 /xm. 

The dashed and dotted lines show the results obtained 
with the CG^ model using recipe 1 and recipe 2, re- 
spectively. Both prescriptions seem to work equally well, 
and improve the SED description at wavelengths longer 
than about 2-3 /im, reducing the discrepancies below ~ 
20%. Note that the larger difference in the Herbig Ae 
SED around 40 /im is not related to scattering (see §3.2). 
Although we do not find any strong reason to prefer recipe 
1 over recipe 2, we note that recipe 2 reproduces more cor- 
rectly the overall reduction of the infrared excess for large 
values of oj (compare the run of with lj with the results 
of Fig.l in DN03). Finally, it is interesting to notice that 
recipe 1 tends to suppress the 10 /im feature with respect 
to the continuum, while this is not the case for recipe 2. 

5. Summary and conclusion 

In this paper we analyze the reliability of two-layer 
models for passive irradiated flaring dusty disks around 
pre-main-sequence stars by comparing their SEDs with 
those computed from more realistic 1-f ID vertical struc- 



100 








Herbig Ae (Fl ) 


Ignore scot 
Recipe 1 


50 




Recipe 2 






/ ■ 





. / i ■■. /\ 




-50 


■ ' 








&; = 0.5 


100 







1 10 100 1000 

A [,iim] 



S 
'q 



T Tauri (F2) 




Ignore scat 
Recipe 1 






Recipe 2 


— j / 






/ ' ^ 






/■' '' 
/.-■/ 







1 10 100 1000 

X [ftm] 

Fig. 13. The factor by which scattering suppresses the 
SED of the disk annulus of test cases Al and A2 for an 
albedo oi ui — 0.5. The solid line is the suppression factor 
from the 1-D vertical structure model. The other lines are 
for the two recipes proposed for the CG+ model. 

ture models. The two-layer model used in this paper is 
called CG+, and is an improved version of the Chiang 
& Goldreich model (see DDNOl). The I+ID model 
involves detailed 1-D vertical radiative transfer which 
treats scattering of the primary stellar photons penetrat- 
ing the surface layers of the disk (see DN03). All the 
models presented in this paper can be downloaded as 
ASCII tables from http : //www . mpa-garching . mpg . de/| 
PUBLICATIONS/DATA/radtrans/cgcompare/. 

We start with a comparison of the SEDs for differ- 
ent disk models where scattering is ignored. In general, 
CG"*" models over-predict the emission at short wave- 
lengths. This is caused by the assumption that the surface 
grains have the temperature of grains in an optically thin 
medium. In many cases one can obtain a better agreement 
by computing the surface temperature with at an attenu- 
ation exp(-0.7) of the stellar flux. At longer wavelengths 
(mid and far-infrared), the CG"*" SEDs sometimes shows 
two maxima, due to the two components of the model: one 
at mid-infrared wavelength (surface layer) and one at far 
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infrared wavelengths (interior) (see, for example, Fig.l^. 
The I+ID model, on the other hand, always produces a 
smooth SED (with the exception of dust features), which 
looks qualitatively somewhat different. However, a closer 
look shows that the agreement between different models 
is often quite good (within about 20%), with some poten- 
tial problems deriving from the use of mean opacities in 
CG+ models. In our set of models, the largest discrepan- 
cies occur for the Herbig Ae star, where they reach ~40% 
around 40 /im. The shape and strength of the 10 micron 
emission feature is in general well predicted by the CG^ 
models, and so is the peak-over-continuum ratio of this 
feature and its intensity profile. We therefore argue that 
also other dust features will be reasonably well reproduced 
by CG+ models. At millimeter wavelengths, both models 
predict similar fluxes (within 20% at most) and intensity 
profiles. 

A comparison of the disk physical structure is also in- 
teresting. Both the pressure scale height Hp{R), which 
depends on the midplane temperature, and the surface 
height Hs{R) predicted by the CG+ model agree reason- 
ably well with the 1-l-lD models. 

Note that we have used in this paper an improved 
version of two-layer code, which does not make use of 
the analythical expressions for flaring angle and temper- 
atures derived in CG97, but derives the disk structure 
self-consistently. A comparison of the CG97 analythical 
model with the results of ID calculations can be found in 
Kraus & Kriigel (in preparation). 

CG+ models do not include scattering of the stellar 
light from dust grains. When compared to the more ac- 
curate 1-l-lD models which include scattering, even for 
quite low values of the albedo or rather peaked-forward 
scattering, we find that the differences in the SED can be 
larger than a factor of two in the near infrared. However, 
as the effect of scattering itself decreases at longer wave- 
lengths, so do the differences between models. We have 
implemented two different simple recipes to our CG"*" code 
to include scattering without changing the basic idea on 
which two-layer models are built. Both recipes seem to 
improve somewhat the agreement of CG+ models with 
the results of the 1-l-lD calculations, although recipe 2 is 
probably a better approximation for very large values of 
uj. (see DN03). 

In summary, the results from this paper show that 
CG+ models compare well with the much more detailed 
1+lD vertical structure models. These two-layer models 
can therefore be used with reasonable confidence to fit ob- 
served SEDs of actual objects, provided that the param- 
eters are not stretched too far. Problems could occur, for 
instance, when the disks are not optically thick enough, 
so that the annulus-by-annulus approach, on which both 
the CG"*" and the 1-l-lD models are based, breaks down. 
In such a case one must resort to a full 2-D or even 3-D 
treatment of the problem. A future study should find out 
under which circumstances such a 2-D/3-D approach is 
unavoidable, or, oppositely, when the annulus-by-annulus 
approach (and thereby the CG"*" approach) is justified. 
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